%q=quad(@myfun,0,1)
%function y = myfun(x) 
    function y1 = myfun(x) % Compute the polynomial.

 % Compute the polynomial.
    %global pp
    global pol_der
	y1 =sqrt(1+(polyval(pol_der,x)).^2);
   % end
    %y1 =polyval(pp,x);%q1=quad(@myfun,0,1);
%y=quad(@myfun,0,x)./q1